Observations of Microwave Continuum Emission from Air Shower Plasmas 
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We investigate a possible new technique for microwave detection of cosmic ray extensive air showers which 
relies on detection of expected continuum radiation in the microwave range, caused by free-electron collisions 
with neutrals in the tenuous plasma left after the passage of the shower. We performed an initial experiment at 
the AWA (Argonne Wakefield Accelerator) laboratory in 2003 and measured broadband microwave emission 
from air ionized via high energy electrons and photons. A follow-up experiment at SLAC (Stanford Linear 
Accelerator Center) in summer of 2004 confirmed the major features of the previous AWA observations with 
better precision. Prompted by these results we built a prototype detector using satellite television technology, 
and have made measurements suggestive of detection of cosmic ray extensive air showers. The method, if 
confirmed by experiments now in progress, could provide a high-duty cycle complement to current nitrogen 
fluorescence observations. 



o 



> 

00 

in 

o 
o 



- 

X 



I. INTRODUCTION 

The origin and nature of the ultra-high energy cosmic rays 
(UHECR) remains one of the enduring mysteries of experi- 
mental particle astrophysics. In spite of well over four decades 
of observations of lO^'' eV UHECR by many different ex- 
periments [T], we still do not have a confirmed astrophysical 
source for these particles, nor do we understand their compo- 
sition in any detail, nor do we know how they propagate from 
their unknown sources to earth f?]. In the last decade new ob- 
servatories such as HiRes and most recently the Auger Obser- 
vatory have much improved the statistics on measurements of 
these particles, but the issues of their origin and propagation 
remain largely open. As the highest energy subatomic parti- 
cles observed in nature, UHECRs must arise from the most en- 
ergetic phenomena in our universe. Their study is thus crucial 
to understanding the nature of acceleration processes that can 
attain energies some seven orders of magnitude higher than is 
currently achievable in the laboratory 1 3, 4, 5,] . 

The primary UHECR spectrum is described by a sim- 
ple power law J{E) oc £-« with a ~ 2.7 for 10^^^ <E < 
2Q19.5 gY Above 10'^-^, the interaction length of cos- 

mic ray nucleons on the cosmic microwave background be- 
comes comparable to intergalactic separation distances, a pro- 
cess first described by Greisen |8], and Zatsepin & Kuzmin |9] 
and now known as the GZK process. It is precisely at and 
above the GZK energies that the measurements of the primary 
UHECR become uncertain due to low statistics, and the shape 
of the spectrum near the endpoint is still a subject of active 
debate. 

Because of the scarcity of particles at these highest en- 
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FIG. 1: "World ultra-high energy cosmic ray and predicted cosmo- 
genic neutrino spectrum as of early 2007, including data from the 
Yakutsk [13], Haverah Park iH/ the Fly's Eye [l6], AGASA [^], 
HiRes [18], and Auger [2dll . collaborations. Data points represent 
differential flux dI{E)/dE, multiplied by E^. Error bars are statisti- 
cal only. GZK neutrino models are from Protheroe & Johnson ^2^1 
and Kalashev et al. [21!]. 



ergies, research into new methods has focused on indirect 
means of observation ll24l Izsll . which makes use of radi- 
ated air-fluorescence emission from the air shower to observe 
it at distances of up to tens of km from the particle axis. 
By observing the longitudinal and transverse development of 
UHECR-induced extensive air showers (EASs) investigators 
are gaining information on the primary composition, which fa- 
vors Ught elements and disfavors a significant electromagnetic 
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(e.g. photon) component. Such studies can also elucidate the 
high energy physics of the early interactions, which occur at 
center-of-mass energies well above that currently probed by 
accelerators [2]. 

The region near the endpoint of the UHECR energy spec- 
trum is shown in summary form in Fig.[Tl where no effort has 
been made to correct the systematic offsets in the flux lev- 
els of the different experiments involved. Above 10^" eV, the 
event rate is of order 1 per km^ per century, producing still 
only a handful of events per year close to this threshold in 
all existing UHECR observatories. As is apparent from the 
current world spectrum, constraints on the high energy tail or 
statistically compelling details of any putative cutoff above 
the current highest energies will still require years of obser- 
vation. The need for much-improved statistics to address the 
primary issues currently under investigation all argue for ex- 
pansions of and improvements on existing methods. Such is- 
sues include the detailed shape of the UHECR energy spec- 
trum |6, 7] (including the presence, or lack thereof, of the 
GZK cutoff [8, 9]), energy-resolved primary particle compo- 
sition [10, |11], a nd source production mechanisms (i.e. ori- 
gins) [I2.1HIII2I. 

In addition, the virtual certainty of the extragalactic origin 
of these particles ensures an associated cosmogenic neutrino 
flux, generated via photohadronic processes throughout the 
universe lfT9ll . Hadronic cosmic rays above ~ lO'^ eV prop- 
agating in the 2.7K cosmic microwave background radiation 
(CMBR) exceed the threshold for resonant A+ particle pro- 
duction through the GZK process, and the rapid decay of these 
unstable secondaries leads to pions and subsequently neutri- 
nos. The mean free path of a 10^" eV proton in the CMBR 
is several Mpc in the current epoch, whereas the neutrinos are 
unattenuated from any cosmic distance. Future observations 
of cosmogenic GZK neutrinos will provide a unique and com- 
plementary view of the UHECR production, propagation, and 
attenuation throughout the universe, motivated by the UHECR 
observations themselves. Figure [1] shows also a band indicat- 
ing the range of models for these cosmogenic neutrinos. The 
uncertainties in the models stem primarily from the details of 
the highest energy part of the UHECR spectrum, as well as the 
epoch of maximum UHECR source evolution, and GZK neu- 
trino observations will thus provide independent constraints 
on the UHECR sources. 



estimates, which may disagree depending on the model used. 
While the optical fluorescence method enables one to observe 
longitudinal as well as transverse shower development, it is 
highly constrained by the fact that it can only work on clear, 
moonless nights. This leads to a net yearly duty cycle of only 
5-10%|241. Furthermore, because the highest energy events 
are observed at increasingly large distances, even small fluc- 
tuations in atmospheric aerosol contamination can have sub- 
stantial effects on energy estimation. 

An air shower dissipates virtually its entire energy bud- 
get through ionization, producing a tenuous plasma with an 
electron temperature of order 10^ K or more. The ioniza- 
tion and subsequent de-excitation of molecular nitrogen in the 
N2^* IN and 2P states leads directly to the optical N2 fluores- 
cence now observed. The hot air shower plasma cools rapidly 
on 1-10 nanosecond time scales, distributing its thermal en- 
ergy through collisions with the neutral molecules, primarily 
N2, which has the largest cross section and number density. 
This rapid cooling process leads to additional excitation of ro- 
tational, vibrational, electronic valence, and other modes of 
kinetic energy distribution among molecules, many of which 
can also lead to rapid de-excitation and subsequent emission. 

In turn, the hot electrons themselves, while producing this 
excitation, can produce their own emission, such as contin- 
uum bremsstrahlung emission, or recombination radiation. 
The fraction of total radiated energy in optical fluorescence, 
compared to the total available energy budget for secondary 
radiation, is very small, leaving much possible radiative en- 
ergy still unaccounted for. The possibilities for observing sec- 
ondary air shower plasma emission other than optical fluores- 
cence have not yet been explored in any detail. 

We report here on investigations of the feasibility of other 
channels for EAS observations. To this end we have under- 
taken several experimental efforts, including two accelerator 
experiments. The promising results from these measurements 
have led us to commission a testbed prototype detector, which 
has helped to establish the methodology that could be used 
to make detailed measurements of EAS microwave molecu- 
lar bremsstrahlung radiation (MBR)|28]. In this report we 
describe the accelerator results and the testbed development 
that has resulted from them, which we have dubbed the Air- 
shower Microwave Bremsstrahlung Experimental Radiometer 
(AMBER). 



1. Motivation for Microwave EAS Detection. 

While there is general agreement among the different ex- 
periments as to the global properties of the UHECR spectrum, 
as Figure [T] shows, there is still significant disagreement and 
uncertainty on absolute flux scales and on some of the fun- 
damental questions of UHECR research. The two primary 
techniques of UHECR observation, ground-based particle ar- 
rays and optical fluorescence detectors both suffer from tan- 
gible limitations. In the case of ground arrays, only a sin- 
gle slice of EAS longitudinal development can be observed. 
This means that determinations of primary particle energy 
and composition require extrapolation via model-dependent 



A. Molecular Bremsstrahlung Radiation. 

MBR in weakly ionized air is created by free electrons ac- 
celerating through collisions with the fields of molecules in 
the ambient medium. EAS ionization is considered "weak" 
since the interactions of free electrons or ions are dominated 
by collisions with neutrals rather than other ions. MBR has 
been classically treated as a thermal process, with the emis- 
sion coming from < 10 eV electrons assumed to be distributed 
with isotropic Maxwellian velocities. By these assumptions 
steady-state MBR emission is expected to be isotropic and un- 
polarized, which strongly differentiates it with the highly di- 
rectional bremsstrahlung from relativistic particles which may 
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FIG. 2: Top: Cross section for momentum transfer between electrons 
and N2, O2, and air over the energy range of interest for molec- 
ular bremsstrahlung production u8l pSI/ . Bottom: The red-dotted 
lines show the function ~E^j^ which must be greater than either 
<5m {weak criterion; blue line) or 2(5 m (strong criterion, black-dash 
line) as a necessary condition for stimulated bremsstrahlung emis- 
sion in an air plasma. See text for details. 



be more familiar to a high energy or cosmic-ray physicist. 
MBR emission shares this property with optical fluorescence 
emission, an important feature which allows for the possibility 
of performing shower calorimetry by mapping the MBR in- 
tensity (and thus the ionization content) along a shower, much 
as optical fluorescence detections maps the nitrogen excita- 
tion along the same shower As indicated in Fig. [TO] MBR 
emission and optical fluorescence both are emitted in all di- 



rections around an EAS, and detectors may therefore "image" 
the air shower glow as a track along the sky to establish two- 
dimensional angular geometry, and use the timing information 
for the pulse arrival to determine the range evolution of the 
shower, giving a complete reconstruction of both geometry 
and particle number evolution. 

The expected isotropic behavior of MBR is also in sharp 
contrast to relativistic radio emission processes such as 
Cherenkov, transition, or synchrotron radiation, all of which 
are beamed and highly polarized. In this respect it is conve- 
nient to think of MBR emission as analogous to "radio fluores- 
cence," whereas beamed relativistic radio emission is closely 
aligned to the particle content of a shower and thus should be 
identified as concordant with the information derived from a 
ground EAS detector array. Furthermore, MBR intensity is 
expected to be proportional to the EAS ionization rate, which 
is known to be itself proportional to A^, the total number of 
charged particles in the shower This therefore leads to a di- 
rect relationship between MBR intensity and shower energy, 
with the degree of proportionality determined by the details 
of local correlations between electron velocities or radiative 
transitions in the tenuous air shower plasma. 

The proportionality will depend on important details which 
require empirical determination, much the same way that os- 
cillator line strengths necessary to understanding optical flu- 
orescence must be determined via laboratory calibration with 
additional corrections for atmospheric conditions. For exam- 
ple, since an EAS produces an initial distribution of ionization 
which is likely to be a power-law rather than a Maxwellian, 
there are corrections for non-thermal effects such as stim- 
ulated emission and other non-equilibrium continuum radio 
emission channels, which may significantly increase the emis- 
sion power over the minimal thermal MBR baseline. In addi- 
tion, the cross sections for both elastic and inelastic collisions 
of electrons with air molecules are complex functions of elec- 
tron energy, yielding strong velocity dependence in the elec- 
tron collision frequency which can also contribute to MBR 
emission coefficients. Such effects are difficult to analytically 
estimate and will be best calibrated in situ as has been done 
for other EAS observation methods. 

To analytically determine the expected minimum flux den- 
sity for MBR, we start with the emissivity from classical 
bremsstrahlung analysis of collisions between electrons and 
neutral molecules [i28i1 : 



ria)(u) = 



(u)C(v. 



,co) 



(1) 



Where co is the microwave radian frequency, u is the elec- 
tron velocity, Ve„ (u) is the velocity-dependent electron-neutral 
collision frequency, and ^(Ven,tt)) is a term that accounts for 
the collisional suppression of radiation which arises from the 
destructive interference of fields from successive collisions 
within the radiation formation zone of each collision, a pro- 
cess also known as plasma dispersion [28]. 

Under the assumption of an isotropic and time-stationary 
velocity distribution, 

^^^-"^^ l + (vju)/co)^ ■ 



4 



For an altitude of 5 km, Ve„ ~ 3 THz at electron energies of 
about 2 eV, near the peak of the collision cross section; for 
room temperature electrons, Ve„ ~ 40 GHz. The correspond- 
ing suppression factors are ^ ~ 5 x 10^^, and C, ~ 0.4, respec- 
tively, showing the wide extremes of values possible under the 
range of electron temperatures that obtain in an air shower 

To preface further discussion below, we note that there are 
several other effects that compete with plasma dispersion and 
will tend to enhance the emissivity, or "suppress the suppres- 
sion." First, stimulated emission, even at a very small level, 
leads to correlations in electron-photon transitions. Such 
emission does not require a full-blown population inversion, 
as in laser processes. Rather, deviations from the ground-state 
Maxwellian distribution can enable low levels of stimulated 
emission. Second, velocity correlations of the electrons due to 
the imposed geometry of the shower tracks and the anisotropic 
distribution of ions can impose some cylindrical symmetry in 
the distributions, in tension with the assumption of uncorre- 
cted isotropy in the electron behavior. Finally, weak plasma 
correlations on the scale of the Debye length can lead to co- 
herent motion of electrons over very small scales, but large 
enough to overcome some fraction of the suppression effects. 

To introduce an ad hoc accounting for the sum of all such 
effects, we impose a coherence factor t, which modifies the ^ 
term: 



density 1^: 



C,-(v,„,co,^) = 



1 + (V,„(U)/C0)2 



(3) 



where ^ > 1 parameterizes the level of excess emission above 
the "suppression floor" determined by ^ in the absence of any 
correlation of either electrons or electron-photon transitions. 
The term t, can then be determined empirically to set the scal- 
ing of the emission, with t, = 1 representing the MBR floor 
level. 

The emissivity r|a) above must be integrated over the distri- 
bution of electron velocities to yield the emission coefficient 
ja (W m^^ (radian/s)^' steradian^*): 



jm= / ria,(u)/(u)u du, 
Jo 



(4) 



where /(u) is the electron distribution function, which is 
Maxwellian in the thermal Umit: 



fW = ( ^^ZITF ) exp 



iTlkT. 



2kT, 



(5) 



for electron temperature Tg. Similar analysis yields the ab- 
sorption coefficient am: 



am = -- 



471 % 

3?co2 



,(u)C(v,co)^^u^t/u 
ou 



(6) 



where the plasma frequency is given by co^ = A^^e^/ (w£,eo) for 
electron number density Ne ■ These coefficients are combined 
to form the source function ~ (l/n^)i(o/ctK) which is then 
integrated along a ray s through plasma column to the observer 
to determine the net intensity per unit radian frequency, or flux 



/m = 



Jo 



e-^dz 



(7) 



where the optical depth T is defined by dz = —a^ds. Note 
that the absorption coefficient is not necessarily positive def- 
inite: under conditions where > 0, e.g., if there is an 
inflection in the electron velocity distribution function, then 
Ua can become negative and stimulated emission will cause 
the radiation to grow with propagation distance. 

The MBR flux density received by a ground detector is es- 
timated by integrating the intensity thus derived over the solid 
angle of the receiving antenna beam to yield Watts per m^ per 
Hz over a given frequency band. The minimum detectable 
change in flux density for a radio antenna and receiver is Ii33ll 



A4 



kBTns 



Ae//\/AfAv' 



(8) 



where Tiyj is the noise temperature of the receiver system, 
Aeff is the effective area of the antenna. At is the receiver sam- 
pling time constant, and Av is the receiver bandwidth. After 
passage of the relativistic shower front, the emission contin- 
ues during the electron thermalization time, f;/,. This time is 
determined by both elastic and inelastic collisions of electrons 
with air molecules 136, 371 , giving f,/, ~ 10 ns for dry air at 1 
atmosphere. The received radiation continues during the time 
the shower remains in the field-of-view of the antenna. At 4 
GHz, for a D=1.8 m dish (giving a beamwidth 2X/D ~ 5°), 
this time can be many microseconds for a distant shower 

Note that there is also a direct analog possible between 
an optical fluorescence detector which used photomultiplier 
"pixels" to image the two-dimensional projection of the opti- 
cal fluorescence intensity along the shower, and a radio dish 
which can use an array of focal plane receiver feeds (each 
of which is effectively a pixel) to image the MBR intensity 
along the same shower. Feed pixel-arrays are not so common 
in radio astronomy because of the success of radio interferom- 
etry, but they are established technology where low-resolution 
bolometric imaging is important (for example, in submillime- 
ter radio astronomy). 

Based on the parameters assumed above, we have numeri- 
cally integrated the flux density for air showers at a distance 
of 10 km, and we find that the MBR floor level of emis- 
sion, including the full suppression term given in equation ?? 
above with ^ = 1, gives an average detection energy threshold 
of order lO'^ eV. However, as we have already noted, esti- 
mates of the MBR emission of UHECR air showers using the 
standard thermal electron formalism here indicate that the en- 
ergy threshold for detectability of the emission may depend 
strongly on the behavior of the modified suppression term t^c, 
which is affected by several classes of non-equilibrium con- 
ditions that can obtain in an air shower plasma. We address 
these in the next section. 
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1. Departures from Equilibrium conditions. 

A summary of the conditions under which departures from 
the MBR floor are expected is given in Table J] Each of these 
conditions may play a role in air shower emission, and we 
discuss each of them in turn. 



TABLE I: Assumptions for the standard MBR derivation compared 
to actual conditions in air showers. 



Assumed Condition, standard 
MBR derivation 


Actual conditions for air 
shower plasma 


Maxwellian (thermal) elec- 
tron speeds 


Non-thermal, cascade power- 
law with high-energy tail 


Isotropic velocity and mo- 
mentum distribution 


Linear ion-trails introduce 
first-order anisotropy 


Time-stationary, in thermal 
equilibrium 


Highly non- stationary, fast- 
transient relaxation 


Collision frequency a weak 
function of electron speed 


N2 cross section a strong 
function of electron speed 



a. Stimulated Bremsstrahlung. As an example of the 
departure from the assumptions regarding the velocity de- 
pendence of the electron collision frequency. Figure |2jtop) 
shows the experimental electron-molecular nitrogen momen- 
tum transfer cross section Om in the energy range of inter- 
est. Most notable is the 2.3 eV resonance due to elastic col- 
lisions that lead to rotational excitation |38]. This resonance 
region is in fact the complement of what is observed in op- 
tical fluorescence-the energy transferred in this region of the 
cross section is released in part through optical fluorescence, 
and the complexity of it is in part mirrored in the structure of 
the optical nitrogen fluorescence transitions. Such highly non- 
linear changes in coUisional parameters with electron energy 
strongly depart from the assumptions used above, and similar 
effects have been found to lead to stimulated emission even in 
highly collisional plasmas under some conditions |28]. 

In fact, stimulated bremsstrahlung emission from gas dis- 
charge plasmas was observed in the 1980's in a series of ex- 
periments 1 29, 31? ]. Necessary, though not sufficient, criteria 
can be stated for stimulated bremsstrahlung emission: 



(9) 



for anisotropic electron distributions whose direction is par- 
allel to the direction of the electric vector of the propagating 
radiation, and 



dE 



(10) 



for isotropic electron distributions. Figure |2jbottom) evalu- 
ates this condition for the molecular nitrogen case shown in 
Fig. 12 top), and it is evident that both of the conditions above 
are strongly satisfied in the neighborhood of the resonance. 
Under such non-equilibrium conditions, an electron popula- 
tion inversion in the ionized region is possible, and this can 



lead potentially to stimulated emission (iS]. Such inverted 
populations have been observed in discharge experiments in 
molecular nitrogen plasmas |34, 35]. 

b. Oxygen attachment. Molecular oxygen has a momen- 
tum transfer cross section significantly lower than that of ni- 
trogen over this energy range, as shown in Fig. |2top), and, 
after weighting for abundance, its effect on the overall mo- 
mentum transfer cross section in air is minimal. Although O2 
does not contribute much to the thermalization of electrons in 
air, it does however play an extremely important role in re- 
moving free electrons once they have thermalized, since the 
three-body attachment cross section to O2 rises steeply at low 
electron energies. In fact the attachment time scale for room 
temperature electrons is comparable to the ^ 10 ns thermal- 
ization time scale for hot electrons in 1 -atmosphere air ll4l]l . 
Once attached to ions, the electrons can no longer contribute 
to the bremsstrahlung continuum radiation. And since the 
initial spatial distribution of the oxygen ions is highly struc- 
tured, and the ions are almost stationary in the short period 
during electron attachment, they impose a rapidly develop- 
ing anisotropy in the removal of electrons from the free dis- 
tribution, creating MBR in the free-bound transition. This ef- 
fect also imposes anisotropy in the velocity distribution of the 
electrons. 

c. Plasma correlations. The Debye length A-o, over 
which an electron is fully shielded from a neighboring ion in 
a plasma, is given by 



A-D = \JtQkTe/{e^-ne) 



(11) 



where n^ is the electron number density. For an E AS plasma at 
an energy of 10 EeV or more, Xo ^ \ —2 cm. Since over this 
distance there may be several hundred ion pairs along each rel- 
ativistic through-going particle track, along with several hun- 
dred to several thousand tracks per square cm in the vicinity 
of the EAS core, electrons do not behave entirely indepen- 
dently but are subject to weak bulk plasma effects at some 
level which will produce phase-space correlations. For our 
case we may class these effects together with the attachment 
effects above; in either case the end result will be parameter- 
ized via equation[3] 



2. Radiative Coherence. 

While the field strength for a single electron is accurately 
described by the MBR theory, the summation of these fields in 
the presence of correlated velocities can significantly alter the 
resulting ensemble field strength. Such alterations, which may 
be produced by intrinsic shower geometry, or electron-photon 
correlations from stimulated emission, or by other plasma ef- 
fects, still may be described via simple vector sums of the field 
strength of each of the radiating particles involved. 

For individual emitters the resultant field strength will grow 
as a phasor sum [42]: 



^ £1 (v) exp {—ik-Xj) 



(12) 
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where Ne is the total number of electrons in the plasma, el (v) 
is the field radiated from a single electron, k is the wave vector 
of the radiation, and x, is is the position of the jth electron 
with respect to the observation point. The total radiated far- 
field power per unit area P/A is given by the magnitude of the 
Poynting flux 

P/A = \SuA = \E?I^ (13) 

where Zq ~ 377 £2 is the impedance of free space. In the limit 
of complete coherence, the phase factors k ■ x, are all unity, 
|£| = A^f. El , and the total coherent power is P^o/, = Pi, 
where Pi is the power radiated from a single electron. Since 
Ne is proportional to shower energy, the coherent power de- 
pends quadratically on the energy of the primary particle. In 
the incoherent limit, the sum of the phase factors corresponds 
to a two-dimensional random walk in the real and imaginary 
components of the resultant field strength, and the total power 
grows as Pi„coh = Ne Pi- 

While in general the partially coherent case requires a de- 
tailed knowledge of the electron phase space distribution func- 
tion, we can get a qualitative sense of the behavior by consid- 
ering a case where the Ne electrons consist of M subgroups 
of ijg electrons each, such that Ne — M/jg. Assume that the /jg 
electrons in each subgroup radiate coherently, but that the sub- 
groups themselves are uncorrected. Thus, while the radiated 
fields from the M subgroups add incoherently, the subgroup 
electrons themselves radiate coherently, and the resulting par- 
tially coherent power is Ppart = 

now quadratic in jig 

rather than Ne- The ratio of the partially coherent power to the 
incoherent power is proportionally 

1^ ^ = ... (14) 

Pincoh Ne Pi 

Similarly the ratio of coherent-to-incoherent power grows as 
Ne- Since the plasma density of ionization electrons in a 
shower scales linearly with shower energy, both the coher- 
ent and partially coherent regimes will yield radiated power 
that grows quadratically with shower energy. In fact, as soon 
as > 10, coherence begins to dominate over the incoherent 
component by an order of magnitude or more. Even modest 
correlations among the shower ionization electrons can thus 
rapidly lead to much larger detected emission than expected. 

We have parameterized these effects using the correction 
term ^ which modifies the collisional decoherence factor ^ as 
described above. In practice empirical data will be required to 
establish the emission constants associated with these factors, 
as is the case for all other emission mechanisms in a real air 
shower. 

B. Accelerator beam tests. 

Motivated by the fact that even the floor level of fully- 
suppressed emission from the MBR process appeared to us 
to be detectable under air shower plasma conditions, we have 
performed two accelerator tests designed to measure the MBR 
in a laboratory air-shower plasma. In these experiments we 



have found good evidence for microwave continuum emis- 
sion with characteristics suggestive of a major departure from 
the standard incoherent MBR emission scenario, not an un- 
expected result given the variety of different non-equilibrium, 
non-thermal, and partially coherent processes that are possi- 
ble. We detail these results here. 



Copper Faraday box AWA INCOBREMS experiment 




anechoic RF absorber 



SLAC T471 experiment 



anechoic RF absorber 




Copper Faraday box 



FIG. 3: Schematic of AWA INCOBREMS (top) and SLAC T471 (bot- 
tom) experiments, which used electron beams to shower in either 
Tungsten or alumina targets to produce ionization inside an anechoic 
Faraday chamber, observed by internal antennas. 




FIG. 4: Views of the exterior and interior of the Faraday anechoic 
chamber used for measurements of microwave continuum emission in 
the INCOBREMS and T471 experiments. The box is approximately 
aim cube in dimensions. 
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FIG. 5: Phase-stable (partially coherent) component for 20 GHz 
emission observed from the 5mm tungsten target, where background 
subtraction of beam-scraping backgrounds was possible. Upper 
pane shows blue as total emission and red as the background due to 
stray ionization from the beam; the lower pane show the background- 
subtracted results. The data is averaged over several thousand beam 
shots. 



1. AWA INCOBREMS. 

In June 2003, the INCOBREMS experiment was performed 
at the Argonne Wakefield Accelerator (AWA). The beam used 
for the experiment consisted of 12 MeV electrons, pulsed in 
~7 nC charge bunches of 1 .2 cm longitudinal thickness (40 ps 
duration) containing ^4 x lO^'' electrons. This gives a typ- 
ical total beam energy of 5 X 10" eV. The beam was col- 
lided into a fixed radiator of tungsten with a thickness of 2 
or 5 mm (tungsten has a radiation length of 3.5 mm), pro- 
ducing a core of photons with energies 5-10 MeV as well as 
some lower energy electrons. Typically 40-90% of the total 
energy was extracted into photons that traversed the cham- 
ber, depending on the radiator. The conversion was necessary 
since the 12 MeV electron energy was well below the criti- 
cal energy in air, and thus inadequate to produce a fully de- 
veloped shower By converting to gamma-rays we avoided a 
large excess negative charge associated with the passage of 
the electron beam through the Faraday chamber The photons 
entered an air-filled ^ 1 m^ copper anechoic Faraday chamber 
which prevented interference from outside electromagnetic ra- 
diation and absorbed transition radiation caused by beam ef- 
fects on the copper C (3.4^.2 GHz), Ku (10.7-11.8 GHz) 
and Ka (20.2-21.2 GHz) band commercial radio receivers 
were mounted on the insides of the chamber to measure sub- 
sequent radio emissions. Figure [3] shows a schematic view 
of the general layout for both this and the subsequent SLAC 
experiments, and Fig. |4] shows several views of the anechoic 



Faraday chamber employed in both cases (here shown with 
the antennas and receivers used for the INCOBREMS experi- 
ment). 

The photon bunches in the AWA experiment typically 
deposited about 1 PeV of ionization energy while passing 
through the Faraday chamber. Since the radiation length of 
electromagnetic particles in air is of order 300 m at sea level, 
the deposited energy is of order 



Echamher ^ - -)(5 X lO" eV) 

300 e 



10^^ eV. 



The mean energy required per ion pair is of order 30 eV, and 
there are thus about 3 x 10^^ ion pairs produced in the cham- 
ber for each beam shot. The distribution of the free electrons 
in the plasma is of course much denser than the equivalent 
PeV cosmic ray air shower Most (^ 80%) of the plasma is 
produced in a central cylindrical region through the chamber, 
with a radius of order 25 cm, and a mean plasma density of 
order 10^ e^ cm^^. 

Based on the expectations of MBR, we expected to observe 
emission that was incoherent, with intensity that scaled lin- 
early with beam current. When our initial observations indi- 
cated that the emission appeared to be scaling coherently, with 
intensity going as the square of beam current, we developed 
analysis methods that attempted to separate the two compo- 
nents, taking advantage of the fact that the phase stability of 
the coherent component allows for it to be subtracted from the 
total emission. (We initially adopted the term "phase stable" 
to describe this emission, since the degree of coherence was 
unknown.) To attempt to separate out the various components, 
we used two thicknesses (2mm and 5mm) of tungsten radia- 
tors to convert the electron beam to bremsstrahlung photons, 
since this provided a different bremsstrahlung energy distri- 
bution and beam emissivity, which could modify the relative 
contributions of the two components. Our measurements were 
unable to separate out any significant incoherent signal com- 
ponent in these data, due primarily to background issues, and 
indicating that the coherent or phase-stable component was at 
least an order-of-magnitude larger. 

In Fig.|5](top), we show results of the AWA measurements 
at 20GHz (Ka-band) using the 5 mm radiator, which gave 
the cleanest background-subtracted results. Partially coherent 
emission was observed about 50 ns after beam passage. While 
the presence of the beam gave clear excess power levels, there 
is also considerable apparent "signal" present when the beam 
was blocked with lead just before our system. We found that 
acceleration and production of the electron beam within the 
of order 15 m length of this accelerator required that a major 
fraction (80% or more in some cases) of the electrons were re- 
moved by upstream collimation (a controlled scraping of the 
beam), but without any way to remove the secondary radia- 
tion (mostly hard bremsstrahlung photons) that this produced. 
This led to a high level of background ionization in both our 
chamber and the surrounding vault, leading to doubts about 
the reliability of the results. This is evident in Fig. |5] where 
the backgrounds with the beam blocked can be seen to at 
times exceed even the apparent signal. We confirmed the pres- 
ence of such backgrounds using external ionization detectors. 
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We also checked carefully whether any portion of these back- 
grounds could be due to radio-frequency interference, and we 
confirmed that this was not the case. 

However, we note that our conclusion regarding these back- 
grounds implies that they are actually stray signal, due to the 
unwanted beam albedo (the beam components that caused 
scattered bremsstrahlung due to impacts with the side-walls 
of the beampipe). Thus it appeared to us that the presence 
of microwave emission from the chamber ionization was un- 
avoidable. To further pursue the investigation with a more 
tightly controlled beam, another experiment was scheduled at 
The Stanford Linear Accelerator Center. 

SLAC T471, July 2004 
-1 1 1 1 1 1 1 1 1 1 1 1 1- 



SLAC T471, July 2004 
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FIG. 7: A plot similar to the previous figure, but now using a cross- 
polarized antenna which was insensitive to radiation polarized with 
the electron beam. The dynamic range of the system was now im- 
proved so that the noise level is determined by thermal noise, and the 
detected microwave emission extends out to 60 ns or more, with an 
exponential decay time constant of about 7 ns. The upper and lower 
dashed red horizontal lines indicate the minimum detectable inten- 
sity, as given by equation\8\ for the single-shot case, and the 100- 
shot average. The diagonal dot-dash lines are the two extreme-case 
estimates for MBR emission: the upper case for no net collisional 
suppression and the lower case for maximal collisional suppression 
of the emission, both for the case where the electron thermalization 
time constant is the source of the 7 ns exponential decay observed. 



FIG. 6: Average microwave emission amplitude from 100 beam shots 
taken near shower-maximum in the 2004 SLAC T471 experiment, 
using a broadband antenna that was polarized along the electron 
beam axis, and was thus sensitive to partially coherent radiation di- 
rectly from the relativistic electron shower as it transited the Faraday 
chamber. A strong initial pulse is seen, with rapid decay, followed by 
a second exponential tail with a longer decay. The noise level is in 
this case determined by the limited dynamic range of the oscilloscope 
used, rather than the thermal noise level. 



2. SLAC T471/E165. 

In the following year, a similar experiment, T471, was per- 
formed at the Stanford Linear Accelerator Center The con- 
figuration of this experiment was largely the same as that of 
INCOBREMS, but additional precautions were taken against 
EMI and beam backgrounds, and verified in lab and beam cal- 
ibration tests. This experiment was coordinated to be operated 
just downstream of the E165 FLASH experiment, which was 
used to do precise calibration of air fluorescence for the HiRes 
collaboration [32]. The SLAC T471/E165 experiments also 
used a precisely controlled, 28 GeV electron beam which was 
collided with a target consisting of 90% AI2O3 and 10% Si02 



to make showers with varying particle number, from to 14 
radiation lengths of material. In T471/E165, the 28.5 GeV 
electron bunches were used directly to create the showers 
with no intermediate conversion to photons via a bremsstrahl- 
ung radiator, as this was unnecessary given the high electron 
energy. Bunches with a typical charge of ^ 2 x 10^ elec- 
trons were used, giving a total shower energy of typically 
6 X 10'^ eV, very similar to those used at AWA. 

Figure |6] shows results from measurements of the emis- 
sion over the 1.5-6 GHz band, using an antenna that was 
co-polarized with the electron shower momentum. Here the 
square of the average signal voltage is plotted vs. the time 
after beam entry into the Faraday chamber. The transit time 
for the chamber is about 3.3 ns for the beam. An initial strong 
impulse is observed at the first causal point in time after beam 
entry. This impulse is found to be highly polarized with the 
plane of polarization aligned with the beam axis and Poynting 
vector, characteristic of transition and radio Cherenkov radi- 
ation. Such emission was anticipated, and is damped almost 
immediately due to the microwave absorber (> 30 dB per re- 
flection even at angles of order 55° from normal incidence) 
that covers the interior of the Faraday chamber (the implied 
average time constant for quasi-exponential decay of reflec- 
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FIG. 8: Top: relative integrated energy in the tail of the microwave 
emission, between 15 to 30 ns, as a function of the depth of the 
shower in radiation lengths. The curve shows a Gaisser-Hillas 
shower profile for comparison, which peaks at about 4. 7 radiation 
lengths at shower maximum. By contrast, the microwave tail emis- 
sion shows some early radiation probably due to the initial electron 
bunch. Near shower maximum, the shower charge overcomes the 
beam charge, and the emission appears to follow the shower pro- 
file near shower maximum, though it falls below the particle number 
profile at late times. 



tions is of order 1 .3 ns for this absorber in this geometry). The 
noise level in this plot is dominated by digitization noise, since 
the sensitivity had to be reduced in order to achieve enough 
dynamic range to see the strong initial impulse. 

In Fig.|7]we plot the same measurements made with a cross- 
polarized antenna, which was therefore insensitive to the rel- 
ativistic shower emission, with a 20 dB cross-polarization re- 
jection factor. In this case the strong initial impulse is not 
prevalent though in fact the leading edge is likely to be slightly 
influenced by the -20dB leakage from the other polarization. 
The exponentially-decaying tail of emission extends out to 
60 ns or more, with noise levels now determined by the ther- 
mal noise level rather than scope noise. Based on beam-out 
vs. beam-in and beam-on vs. beam-blocked measurements, 
no beam-related background (either ionization or EMI) was 
present, and thus no background subtraction was necessary. 
Several curves are also plotted with the data. The horizon- 
tal lines indicate the thermal noise level for single shots and 
for the average of the 100 beam shots used here, based on 
equation [8] above. The diagonal dot-dash curves are model 
predictions based on equation|2]above, calculated for the two 
extreme cases of the collisional term from equation[3j one 
for the case of no suppression (^c = 1) and the other for full 
collisional suppression = 1). It is evident that, if MBR is 
responsible for this emission, the collisional suppression is al- 
most completely offset by the partial coherence. 



In Fig. IS top) we plot the behavior of the integrated mi- 
crowave energy in the 15-30 ns window as a function of 
shower depth in radiation lengths. The emission from the 
direct beam, which adds noticeably to the shower emission 
up to about 4 radiation lengths, has been subtracted here in 
proportion to the depth in radiation lengths, so that the con- 
tribution from the shower emission alone can be compared to 
expectations. The upper (long-dashed) curve shown is a stan- 
dard Gaisser-Hillas profile |40], peaking at about 4.7 radiation 
lengths for these showers. The lower (short-dash) curve is a 
Gaisser-Hillas profile, but now scaled as the square of the par- 
ticle number in the shower. It is evident that the shower emis- 
sion scales roughly with the particle number in the shower, 
but appears to drop below the standard Gaisser-Hillas profile 
at large shower depths, although not enough to warrant scal- 
ing that is quadratic in particle number. This behavior pro- 
vides evidence that the process for the emission is relatively 
insensitive to the plasma density. At larger shower depths in 
particular the plasma density decreases by 1-2 orders of mag- 
nitude with only factors of 2-3 apparent drop in the relative 
microwave emission compared to expectation based on the 
standard Gaisser-Hillas profile. 

Fig.|8jbottom) shows that the fitted time constant of the de- 
cay of the emission power is roughly constant with shower 
depth at about 7 ns, with some indication that it may be in- 
creasing for large shower depths. The near constancy of this 
parameter indicates that the underlying physical process that 
removes the radiating electrons from the emitting populations 
is nearly completely insensitive to plasma density. 

The radiation observed in T471 is also partially coherent. 
This is shown in Fig.|9] which plots the integrated microwave 
power from 15-30 ns after the main pulse vs. beam energy 
as measured by an external transition radiation detector. The 
quadratic correlation here indicates that the partially coher- 
ent portion of the emission dominates completely over inco- 
herent emission. The coherent emission appears to be sev- 
eral thousand times the expected incoherent emission level, 
implying that subgroups of /j^ ~ 10^^^^ electrons are radi- 
ating quasi-coherently within their subgroup, using the nota- 
tion of equation [14] above. Given that the showers used in 
T471 created 3 x lO'^ ionization electrons within the Faraday 
chamber used, the net correlation of 10^ is still extremely 
small, and it is evident that this level of partial coherence is 
very far removed from the full-spatial coherence that obtains 
in coherent Cherenkov or transition radiation. 

We note that the Debye length (equation [TT] above) for the 
T47 1 plasma is initially of order 2 mm when Te ~ 10^^ K, and 
within one Debye radius there are of order 10^ free electrons 
initially. Thus a weak correlation of ~ 0.01% within a Debye 
radius appears to be all that is required to create the observed 
partial coherence effects. This analysis does not account for 
the rapid evolution of the Debye length as the electrons cool 
however. At Tg = IQ^ K, close to ambient, Xo ~ 0.3 mm, 
and the Debye volume then contains of order 10^ electrons, 
still requiring only a 1 % correlation coefficient. However, any 
prediction using a simple plasma-correlation model requires 
understanding of the dynamics of the cooling event before a 
self-consistent picture can emerge. 
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FIG. 9: Plot of relative microwave energy in the tail of the observed 
air plasma emission vs. the relative microwave energy observed in 
an external transition-radiation monitor of the beam current, which 
is directly proportional to beam energy. The observed microwave 
power follows closely a quadratic rise with beam energy, character- 
istic of coherent radiation. 



3. Scaling to air showers. 

Under the assumption that standard radiation scaling laws 
obtain, we can make an estimate of the threshold for de- 
tectability of the emission seen in Fig. |2l if we scale it to air 
shower observation distances and a realistic detection system. 

To do this, we take the flux density as estimated from the 
data in Fig.|2l using the weighted average T471 antenna effec- 
tive area of 0.05 m^. Using this flux density and the equivalent 
shower energy of 3.4 x lO'^ eV, we scale to an equivalent air 
shower at a distance of 10 km. The scaling corrects for the 
lower electron density expected for a typical 5 km air shower 
altitude. We also assume an integration time (several hundred 
ns) based on angular transit times for showers at roughly this 
distance, with the peak flux density determined by the emis- 
sion over an interval comparable to the thermalization time. 
We consider both linear and quadratic scaling of the emission 
with electron density, and with regard to the shower geom- 
etry, as long as the transverse diameter of plasma column is 
contained within the antenna beam, we assume there is a di- 
rect scaling from the T471 observations to air shower obser- 
vations. 

The results of this analysis indicate that, if the partially co- 
herent emission observed in T471 scales only linearly with 
shower energy (as might be expected in the pure "incoher- 
ent" case), then the threshold of the AMBER system at 10 km 
is of order 8 x 10^^ eV in shower energy. If the scaling is 
quadratic with shower energy, as the data suggests, the thresh- 
old is lower, of order 1.6 x 10^** eV. 



We can also estimate the maximum distance to which a 
shower at the GZK threshold energy of 3 X 10''^ eV could be 
observed under these same conditions: for the linear-scaling 
case, the distance is of order 20 km; for the quadratic case it 
is much larger, of order 200 km, but of course in this case 
earth-curvature and atmospheric attenuation would also re- 
quire consideration. In either case, the current emission pa- 
rameters strongly warrant further investigation of the potential 
for development of MBR detection of air showers. 

As noted in a previous section, an air shower plasma can 
also be characterized in terms of its Debye length. For a 
lO'^ eV shower, with an initial electron density of order lO'" 
e^ m^-' within a few m of its core, Xo — l cm for Te = \(f K 
during the early period of the electron cooling, and A-d ~ 1 cm 
once the electrons have cooled close to ambient levels. The 
electron number within a Debye radius similarly evolves from 
~ 10^ to ~ lO'* over the same cooling period, several tens 
of ns. These values leave open the possibility of correlated 
electron behavior comparable to those seen in our T471 ex- 
periment, if the correlations are related to plasma density pa- 
rameters. 



4. Beyond detection: Shower Calorimetry with an AMBER array. 

The importance of MBR detection of EAS rests in the po- 
tential that it will yield the observational advantages compara- 
ble to those of optical fluorescence without the shortcomings 
associated with weather and limited duty cycle. By observing 
MBR, one is observing an EAS from the same perspective as 
with optical fluorescence, via energy-loss processes that are 
intimately related to the excitation of molecular nitrogen that 
leads to air fluorescence. However, observations can occur 24 
hours per day, and at the microwave bands of interest there 
is virtually no attenuation due to atmospheric contamination 
from aerosols or clouds. Even heavy rain leads to attenuation 
of< 1 dB above elevation angles of 30°at C-band (4-6 GHz), a 
20% effect. Initially, while the MBR technique is being cross- 
calibrated with respect to an optical fluorescence and ground 
array, this immunity to atmospheric effects can yield imme- 
diate benefits in helping to extrapolate the energy scale for 
distant events, where optical fluorescence is most affected by 
aerosols and other atmospheric uncertainties. 

Commercially designed microwave reception equipment 
can be easily weatherproofed, and future arrays would most 
likely be able to employ off-the-shelf satellite television com- 
ponents, taking advantage of the tremendous economy of 
scale in wireless and satellite television technology. Follow- 
ing validation of the technique in coincidence with an exist- 
ing EAS installation, MBR detectors could be potentially de- 
ployed as standalone UHECR observatories. Critical to the 
success of such an observatory is the ability of MBR to do 
precision shower calorimetry. There appears to be good ini- 
tial evidence from the T471 experiment that such calorimetry 
can be done with precision that is comparable to current tech- 
niques. 

Referring to Figure |9l we stress that individual single- 
shot measurements of the integrated microwave energy in the 
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TABLE IL Detectability of air showers assuming that the T471 shower emission shown in Fig. |7] can be scaled to an air shower at 10 km 
distance, reaching maximum at 5 km altitude, and detected with a 1.8 m diameter dish at 4 GHz (C-band) using standard satellite dish 
technology. Results are estimated both for linear and quadratic scaling with shower energy. 



Inputs and derived parameters 


expression 


value 


units 


Observation parameters 




peak intensity observed, T471, single pol 




1.00E-06 


W/m'^2 




time constant 


T 


1.00E-08 


s 




reference bandwidtli 


Afo = A(o/2n 


2.50E+09 


Hz 




reference system temperature 


Tsys, ref 


4.00E+02 


K 




reference flux density 


IfD = lo/ Afo 


4.00E-16 


W/m^/Hz 


Reference shower parameters 




electron energy 


e 


2.80E+10 


eV 




beam current 


jb 


1.20E+07 


e-/pulse 




observed reference shower length 


Lo 


0.65 


m 




distance to reference shower axis 


d 


0.5 


m 




total reference shower equivalent energy 


Eo = £jb 


3.36E+17 


eV 


Parameters for estimation of trial air shower signal 




Trial distance to air shower 


R 


1.00E+04 


m 




mean altitude of trial shower maximum 


a 


5.00E+00 


km 




Trial dish diameter 


D 


1.8 


m 




feed efficiency 


Tl 


0.8 






dish effective area 


A,ff=ii7i(D/2f 


2.036 


m^ 




center frequency of observation 


fc 


3.80E+09 


Hz 




receiver bandwidth 


Af 


5.00E+08 


Hz 




shower length over 1 e-folding decay time 


L^ = C T 


3.00E+00 


m 




Ratio of electron density at shower altitude to sea level 
electron density in reference shower 


p 


6.00E-01 






Ratio of trial shower length per e-folding time to 
observed reference shower length 


r = L,/Lo 


4.62 






Expected signal flux density 


If.exp=lK)rp(R/d)-' 


2.77E-24 


W/m^/Hz 


Parameters for estimate of minimum detectable flux density 




feed+dish beamwidth on sky 


Ae= 1.22 c/(fD) 


5.35E-02 


rad 




observed average shower length 


La, = RAG 


9.90E+02 


m 




average duration of emission in field of view = 
integration time 


At = Lav / C 


3.27E-06 


sec 




System temperature 


T 

sys 


8.00E+01 


K 




Boltzmann's constant 


k 


1.38E-23 


W/K 




minimum detectable flux density 


AI = kT,y,/[Aeff(AtAf)"'] 


1.34E-23 


W/m^/Hz 


SNR and predicted threshold energy 




signal to noise ratio for reference shower at R 


SNRo = lf,exp/AI 


2.06E-01 






Threshold shower energy required for 5ct, quadratic 
power dependence with shower energy 


E,h,esh=(5/SNRo)'"Eo 


1.65E+18 


eV 




Threshold shower energy required for 5a, linear power 
dependence with shower energy 


Eth,.e.h = (5/SNRo)Eo 


8.14E+18 


eV 




Maximum distance a 3 x 10^^ eV shower could be 
observed, assuming linear scaling 


Rmax = (3xlO" / Eihresh R 


1.92E+04 


m 




Maximum distance a 3 x 10 eV shower could be 
observed, assuming quadratic scaling 


Rraax=(3xlO"/E„„esh)R 


1.81E+05 


m 
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Particle shower impact with ground: 

- Direct Detection of shower 'slice' by ground array 

- Indirect detection of integrated profile via beamed radio synchrotron 



Developing air shower: 

- Indirect detection of profile of ionization density by: 
a) Nitrogen fluorescence (optical) 




AlVIBER array telescope array 



FIG. 10: Conceptual sketch of how microwave bremsstrahlung detection relates to other methods of ultra-high energy cosmic ray air shower 
detection. 



shower can be used to determine the beam energy in the exper- 
iment to about 2-3% precision, once the overall energy scale 
is set (in this case by an external beam current monitor). This 
level of precision is set entirely by the instantaneous signal-to- 
noise of the data, and is not systematics-limited, as evidenced 
by the precise quadratic scaling observed. In the T471 ex- 
periment, this scaling was observed over more than a decade 
(not shown in Fig.|9l) of energy. In addition, the same scaling 
was seen at many different relative plasma densities (created 
by sampling the shower at different depths of development) in 
the experiment, indicating that it is not dependent on shower 
age. Such results give us good confidence that, with suffi- 
cient attention to careful detector and system design, and ade- 
quate calibration, an MBR observatory could provide shower 
calorimetry which was comparable to that of existing tech- 
niques. 



EAS observatories have also demonstrated capabilities for 
neutrino observations, but will require substantial increases in 
their apertures before such techniques can become practical 
in elucidating the GZK neutrino spectrum. Neutrino-induced 
showers are also likely to be highly-inclined relative to typical 
proton showers, and thus become problematic for observation 
with ground arrays, which suffer from more severe system- 
atics at high zenith angles. In contrast fluorescence methods 
(and possibly the MBR methods we describe here) can read- 
ily observe such showers, since the geometry is no less fa- 
vorable for horizontal than for vertical showers. Thus MBR 
observations may help to greatly expand the neutrino aper- 
tures of air shower observatories, by extending the duty cycle 
for "quasi-fluorescence" observations, perhaps by an order of 
magnitude or more. Figure[TO]gives a schematic view of how 
such methodology relates to other implementations of ultra- 
high enery cosmic ray air shower detection. 



C. The AMBER system. 

Following the indications of stronger-than-expected emis- 
sion from the two accelerator experiments detailed above, we 
have moved ahead to develop a prototype of a system that 
could be used to search for detectable microwave emission 
from actual air showers. This system is built around a custom 
compact-PCl digitizer and data acquisition system, which we 
designate the Radio Bremsstrahlung Impulse Detector (Ra- 
BID). We have chosen the components and size of the pro- 
totype system such that it can be duplicated at low cost with 
mostly commercial parts. The proposed system, incorporat- 
ing the RaBlD prototype, is designated the AMBER for Air- 
shower Microwave Bremsstrahlung Experimental Radiome- 
ter AMBER is currently operating on the rooftop of Watan- 
abe Hall at the University of Hawaii at Manoa (UHM) in Hon- 
olulu, Hawaii, pictured in Fig.fTTI 

In its current configuration, the AMBER unit consists of 
a dual-band (C and Ku), dual-polarization feed horn array at 
the prime focus of a 1.8 m off-axis parabolic dish. The ar- 
ray is in a diamond-shaped configuration where each feed is 
~ 5.2° from its nearest neighbor. Each feed produces four 
channels of signal which are amplified and down-converted in 
Low Noise Blocks (LNB) and then conveyed to the RaBlD 
DAQ via RGll coaxial cable, as shown in Figure [TT] The 
RaBlD DAQ consists of a pair of RaBlD cards located inside 
a compact PCI (cPCl) crate, along with a cPCl CPU for data 
collection and logging. At the RaBlD card input, the down- 
converted LNB outputs are measured with RF power moni- 
tor (MAX4003) chips, which provide output proportional to 
the received RF power, with approximately 70 ns integration 
time. This power level is sampled with a 32MSa/s ADC and 
processed inside a Field Programmable Gate Array. 

These digitized samples are processed in 3 parallel paths: 
(1) all C-band samples are logged into a hardware histogram- 
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RaBID cPCI card (x2) 



compact PCI (cPCI) Data Acquisition crate 

FIG. 1 1 : Top: Prototype AMBER telescope and feed array on the roof 
of the physics building at UHM. Bottom: AMBER detector readout 
chain. The feed horn signals are amplified and down-converted in a 
Low Noise Block and then transmitted to a pair of RaBID cards for 
processing. See text for details. 



mer, which allows optimum threshold-riding with varying 
background (2) a trigger threshold is set based upon the his- 
togram values; and (3) a circular buffer holds the samples in 
time sequence to be read out into the CPU upon detection of 
a trigger condition. In order to avoid biases in the trigger- 
ing, each feed horn channel (of 12 total) is triggered sepa- 
rately, at minimum possible threshold, and the trigger times 
(corresponding to different transit times across the array field 
of view) are analyzed in the stored data. All sample times are 
recorded with respect to a common clock, which is synchro- 
nized to GPS via Network Time Protocol. An external trigger 
port (not shown) is available for forcing readout when observ- 
ing in conjunction with another detector. 



Pseudoronge 1.5 km 




01234 01234 
time, microseconds time, microseconds 



FIG. 12: Example of an event recorded recently with the prototype 
RaBID system in Hawaii. The two events occurred sequentially in the 
two feeds noted, and triggered both polarizations at both feed 4 (top) 
and feed I (bottom), creating a downgoing event trigger. Left,: Raw 
data, showing peaks in both H and V polarizations for the C-band 
feeds. Right: Signal after co-adding the polarizations and applying 
Wiener filtering to remove the high frequency thermal noise fluctu- 
ations. The pseudorange is based on the 680 ns dual-feed crossing 
time. 



1. AMBER Results. 

Since initial commissioning of the AMBER system began 
in mid-2005, we have accumulated about 8 months of data 
under stable operating conditions, most of which has been an- 
alyzed to search for EAS-like events. Because AMBER lacks 
a ground-truth EAS array to validate any observed signals, 
any candidates that are observed remain only putative at best. 
However, we may test a sample of such candidates for simi- 
larity to expectations from our simulations, and we have done 
this for a large data sample taken through the spring of this 
year, with results that support the potential for EAS measure- 
ments by an AMBER array. 

Data analysis for AMBER events involves several steps 
which significantly improve the signal-to-noise ratio of the 
raw data. First, because the MBR signal is by nature unpo- 
larized, we can combine the recorded power in the two inde- 
pendent polarizations, improving the SNR by a factor of order 
\/2. Also, the signal arrives over many sequential 70 ns time 
bins, whereas a significant fraction of the noise is broadband 
and largely uncorrected between successive time bins. These 
statements are equivalent to saying that the spectral bandwidth 
of the signal is much less than that of the noise, and under such 
conditions we may apply Wiener filtering (also known as op- 
timal filtering) to remove the out-of-band noise component, 
and properly weight the in-band noise components. 
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FIG. 13: Examples of other events recorded by the prototype AMBER system, meeting the criteria for EAS candidates. (In these data, the 
convention is reversed compared to Fig. \12\ with the lower curve for each event originating from a feed which scans a higher elevation in the 
sky.) 
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FIG. 14: Upper: Pseudorange distributions of simulated events for 
an AMBER array with a 1 EeV energy threshold at 10 km distance 
(upper curve) and a 10 EeV threshold at 10 km (lower curve). Lower: 
Pseudorange distribution of 10 candidate events measured in recent 
AMBER data taken over several months. 



Examples of the effects of this analysis are shown for a can- 
didate event in Fig. [12] On the left side the raw event that 



triggered the system at C-band (4 GHz) is shown, with the 
upper (earlier) signals from feed 4 and the later signals from 
feed 1 (lower left), indicating a downgoing event. On the 
right-side panes, the signals are shown for the combination 
of polarization co-adding and Wiener-filtering, with a marked 
improvement in overall SNR and resulting timing. Although 
the actual range to the event cannot be determined directly, 
we calculate a pseudorange based on the assumption that the 
feed-crossing signals were moving at the speed of light over 
the known angle between feeds. These pseudorange values 
can then be compared to simulations for actual EAS events. 

Over the several months' observation period where the data 
has the highest quality, we have selected a sample of candi- 
date events based on criteria derived from EAS expectations. 
Additional examples of such candidates are shown in Fig.[T3] 
Here the projected elevations and azimuths for each feed are 
shown in each event pair. We find that downgoing events pre- 
dominate in our current sample. This is expected from a true 
EAS sample, but without an independent air-shower tag for 
any given event, we cannot yet reject the possibility of anthro- 
pogenic origin. 

We can however compare the derived pseudorange distri- 
butions for both simulations and actual data to determine if 
the candidates observed in our event sample are drawn from a 
distribution that is consistent with what is expected from ac- 
tual EAS events. To do this we have developed a Monte Carlo 
simulation code from which we can extract the pseudorange 
value for events with various detectability thresholds and an 
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energy spectrum consistent with the known UHECR energy 
spectrum. 

Resuhs from this analysis are shown in Fig. [14] The up- 
per pane shows the simulations for two energy threshold at a 
distance of 10 km, and the lower pane shows the results for 
a sample of current candidate events. While these candidates 
cannot be proven to be EAS events without independent evi- 
dence from an air shower array, they do appear at least consis- 
tent to first order with the expectations from EAS events, and 
they demonstrate that an AMBER array has the basic detector 
characteristics to make measurements that are necessary to es- 
tablish MBR observations as a viable EAS detection method- 
ology. 



In summary, we have proceeded as far as possible with 
MBR studies in the absence of coincident EAS ground-truth 
array measurements. Efforts are now underway to deploy an 
AMBER testbed array within the Auger Observatory in Ar- 
gentina. 
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